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Abstract 

Cephinae is traditionally divided into three tribes and about 24 genera based on morphology and host 
utilization. There has been no study testing the monophyly of taxa under a strict phylogenetic criterion. 
A molecular phylogeny of Cephinae based on a total of 68 sequences of mtDNA COI gene, representing 
seven genera of Cephinae, is reconstructed to test the traditional limits and relationships of taxa. Mono- 
phyly of the traditional tribes is not supported. Monophyly of the genera are largely supported except for 
Pachycephus. A few host shift events are suggested based on phylogenetic relationships among taxa. These 
results indicate that a more robust phylogeny is required for a more plausible conclusion. We also report 
two species of Cephus for the first time from Turkey. 
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Introduction 

The Cephidae is a small family of Hymenoptera with a thin integument, usually black 
or dark colored and commonly with narrow yellow bands on the abdomen. It com- 
prises approximately 1 60 species in three subfamilies and about 24 genera and is pri- 
marily Holarctic (Benson 1935, 1946; Smith and Schmidt 2009; Taeger et al. 2010). 
Two of the subfamilies, Athetocephinae and Australcephinae are represented by only 
four species and are restricted to Madagascar and the Australian Region (Benson 1935, 
1946; Smith and Shinohara 2002a; Smith and Schmidt 2009). The majority of species 
are included in the Holarctic subfamily Cephinae. Although several faunal treatments 
(Ries 1937; Middlekauff 1969; Zhelochovtsev and Zinovjev 1988; Calmasur and Oz- 
bek 2010; Korkmaz et al. 2010) and a single world review (Muche 1981) have been 
published, and a number of cephid species have been described in recent years (Smith 
and Solomon 1989; Smith 1997; Smith and Shinohara 2002a,b; Smith and Scruff 
2005; Wei 2007; Smith and Schmidt 2009; Wei and Smith 2010), their phylogenetic 
relationships have not been investigated. 

Cephidae can be easily identified since they are morphologically intermediate be- 
tween the hymenopteran suborders Symphyta and Aculeata. Because of several apocri- 
tan-like characters, such as a weak constriction between the first and second abdominal 
segments, the lack of cenchri and the rough area on fore wing, and the form of male 
genitalia, they were once considered as a likely sister group of Apocrita (Konigsmann, 
1977). However, considerable evidence from both morphological and molecular data 
strongly support a sister group relationship between Orussidae and Apocrita, and the 
Cephoidea, containing the only family Cephidae, appears as to share a last common 
ancestor with a lineage leading to the Siricoidea (e.g., Rasnitsyn 1980,1988; Basibuyuk 
and Quicke 1995; Vilhelmsen 1997, 2001; Ronquist et al. 1999; Schulmeister et al. 
2002; Sharkey 2007). 

Ries (1937) suggested that Janus is the most primitive genus of Cephinae based 
on its filiform and many segmented antennae and tarsal claws. Benson (1946) di- 
vided the Cephinae into three tribes, Cephini, Hartigiini and Pachycephini based 
both on morphology and host utilization. The larvae of Cephini bore in the stems 
of Poaceae, those of Pachycephini live in the stems of Papaveraceae and Poaceae, 
and those of Hartigiini bore in the twigs of Rosaceae or other arborescent plant 
families (Benson 1946; Middlekauff 1969). Numerous morphological characters 
and color patterns that traditionally have been used for separation of the taxa with- 
in the family are claimed to be either variable or display phenotypic plasticity (Ries 
1937; Benson 1946, 1968; Korkmaz et al. 2010). Current classification is mainly 
based on morphology and host usage and therefore necessitates a close examination 
under the phylogenetic approaches. 

Phylogenetic studies of taxa that exhibit adaptive phenotypic variation provide 
valuable insights into the evolutionary forces driving the origins of diversification 
(Zhang et al. 2008). Research on phytophagous insects has confirmed that adapta- 
tion and specialization to different plant species are central to generating diversifica- 
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tion at all hierarchical levels (Mopper and Strauss 1998; Berlocher and Feder 2002; 
Funk et al. 2002; Nosil et al. 2002; Eubanks et al. 2003; Lozier et al. 2007). How- 
ever, host specialization in the Cephinae, as observed in most phytophagous insect 
groups, might have led to evolutionary shift between higher taxonomic groups. 
Host shift probably has taken place many times at different periods, and therefore 
the classification based on host use may not reflect the true phylogenetic relation- 
ship within the Cephinae. 

Here, we selected the mitochondrial cytochrome oxidase subunit I (COI) gene to 
reconstruct the phylogenetic relationships of the Cephinae and identify systematic po- 
sition of its tribes and genera by applying phylogenetic inference methods. The selected 
COI gene region is informative for estimating relationships at both intra- and inter- 
species level due to possession of both completely conserved and variable regions and 
having a heterogeneous evolutionary rate across the gene (Lunt et al. 1996; Dowton 
and Austin 1997; Caterino et al. 2000; Roe and Sperling 2007; Bacci et al. 2009) This 
region is also utilized as DNA-based bio-identification system for animals at the global 
level (Hebert et al. 2003). Cephus parvus (Dovnar-Zapolskij, 1931) and C. runcator 
Konow, 1896, are recorded for the first time from Turkey. 



Material and methods 

Sixty-eight specimens representing three tribes and seven genera of the subfamily 
Cephinae were collected from localities presented in Table 1. All specimens are de- 
posited in the Entomological Collection of Cumhuriyet University, Sivas (ECCUS). 
A specimen of Arge sp. (Argidae) was included as an outgroup in the phylogenetic 
analyses. Several keys were used to identify the specimens (Benson 1946, 1951, 1968; 
Muche 1981; Zhelochovtsev and Zinovjev 1988). The taxa names, the voucher speci- 
mens, and GenBank accession numbers are presented in Table 1. 



DNA extraction, amplification, and sequencing 

Alcohol-preserved specimens were allowed to dry on filter paper, and DNA was ex- 
tracted from left legs of the specimens using the High Pure PCR Template Prepara- 
tion Kit (Roche Diagnostics, Mannheim, Germany) following the protocol for DNA 
isolation from mammalian tissue. Each DNA sample was dissolved in 200 ul elution 
buffer and stored at -20°C. The partial mitochondrial COI gene (750 bp) was ampli- 
fied by using the conserved COI primers with the following sequence: COI-sl859, 
5' - GGAACIGGATGAACWGTTTAYCCICC - 3' and COI-a2590 5' - GCTC- 
CTATTGATARWACATARTGRAAATG - 3' (Simon et al. 1994). PCR reactions 
were conducted with 10 ul of extracted DNA in 50 ul reaction mixture. Amplification 
conditions were as follows: denaturation for 5 min at 94°C, followed by 37 cycles of 
denaturation at 94°C for 30 s, annealing at 59°C for 45 s, extension at 72°C for 30 s 
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Table I . List of taxa and voucher specimens used for sequencing. 



Genus 


Species 


Voucher no. 


GenBank ac- 
cession no. 


Location 


Col. date 


Arge 


sp. 


ECCUS 201 


JF901916 






Calameuta 




Jilijormis (Eversmann, 1 847) 


VCC\ IS 7 1 0 


ipqni 84Q 


trel 


1 9 n4 9nnQ 




Jilijormis 




iFQ.ni ssn 




n4 nfi 9nnQ 




haemorrhoidalis (Fabricius, 1781) 


FiTl IS 71 ? 


ipqni 8^9 


fvutanya 


9n ns 9nnQ 




haemorrhoidalis 


FPPT IS 913 


ipqni 8S3 


ispai La 


1 7 n^ 9nnQ 

1 / .KJ J .ZUU J 




haemorrhoidalis 


ECCUS 214 


JF901855 


Kocaeli 


04.05.2010 




haemorrhoidalis 


ECCUS 215 


JF901856 


Kocaeli 


04.05.2010 




haemorrhoidalis 


ECCUS 216 


JF901857 


Bayburt 


05.06.2010 




haemorrhoidalis 


ECCUS 217 


JF901858 


Usak 


19.05.2009 




haemorrhoidalis 


ECCUS 218 


JF901859 


Isparta 


17.05.2009 




idolon (Rossi, 1794) 


ECCUS 219 


JF901851 


Konya 


17.05.2009 




pallipesQQug, 1803) 


ECCUS 220 


JF901854 


Sivas 


13.05.2010 




pallipes 


ECCUS 221 


JF901860 


Hakkari 


11.06.2003 




pygmaea (Poda, 1761) 


ECCUS 222 


JF901848 


Hatay 


09.04.2009 




sp. 


ECCUS 223 


JF901861 


Sivas 


17.06.2007 


Cephus 




brachycercus 'Ihomson, 1871 


FPPT IS 730 


1FQ.D1 871 


Istanbul 


n8 ns 9ni n 

UO.U J.ZU1U 




brachycercus 


ECCUS 231 


JF901872 


Sivas 


10.05.2010 




fumipennis Eversmann, 1 847 


ECCUS 232 


JF901873 


Ardahan 


07.06.2010 




nigrinus Thomson, 1871 


ECCUS 233 


JF901874 


Istanbul 


08.05.2010 




parvus (Dovnar-Zapolskij, 1931) 


ECCUS 234 


JF901875 


Sivas 


17.05.2010 




parvus 


ECCUS 235 


JF901876 


Sivas 


26.05.2010 




pulcher Tischbein, 1852 


ECCUS 236 


JF901877 


Erzurum 


06.06.2010 




pygmeus (Linne, 1 767) 


ECCUS 237 


JF901911 


Denizli 


18.05.2009 




pygmeus 


ECCUS 238 


JF901912 


Hatay 


09.04.2009 




pygmeus 


ECCUS 239 


JF901913 


Hatay 


09.04.2009 




pygmeus 


ECCUS 240 


JF901914 


Bayburt 


07.06.2008 




pygmeus 


ECCUS 241 


JF901915 


Bayburt 


07.06.2008 




rjabovi Dovnar-Zapolskij, 1926 


ECCUS 242 


JF901878 


Kmkkale 


20.06.2009 




rjabovi 


ECCUS 243 


JF901879 


Kirikkale 


20.06.2009 




runcator Konow, 1896 


ECCUS 244 


JF901880 


Edirne 


07.05.2010 




runcator 


ECCUS 245 


JF901881 


Edirne 


07.05.2010 




sareptanus Dovnar-Zapolskij, 1928 


ECCUS 246 


JF901882 


Erzurum 


06.06.2010 




sareptanus 


ECCUS 247 


JF901883 


Erzurum 


06.06.2010 




sp. 


ECCUS 248 


JF901884 


Bilecik 


05.05.2010 




sp. 


ECCUS 249 


JF901885 


Bilecik 


05.05.2010 




sp. 


ECCUS 250 


JF901886 


Canakkale 


06.05.2010 




sp. 


ECCUS 251 


JF901887 


Amasya 


02.05.2010 




sp. 


ECCUS 252 


JF901888 


Amasya 


02.05.2010 




sp. 


ECCUS 253 


JF901889 


Tekirdag 


08.05.2010 




sp. 


ECCUS 254 


JF901890 


Sivas 


18.05.2010 




sp. 


ECCUS 255 


JF901891 


Erzurum 


06.06.2010 
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Genus 


Species 


Voucher no. 


GenBank ac- 
cession no. 


Location 


Col. date 




sp. 


ECCUS 256 


JF901892 


Kars 


07.06.2010 




sp. 


ECCUS 257 


JF901893 


Kars 


07.06.2010 




sp. 


ECCUS 258 


JF901894 


Bolu 


04.05.2010 


Tracheitis 




iudaicus (Konow, 1 907) 


ECCUS 260 


JF901865 


Bayburt 


05.06.2010 




iudaicus 


ECCUS 261 


JF901866 


Bayburt 


05.06.2010 




libanensis (Andre, 1881) 


ECCUS 262 


JF901867 


Icel 


13.04.2009 




libanensis 


ECCUS 263 


JF901868 


Icel 


13.04.2009 




sp. 


ECCUS 264 


JF901862 


Sivas 


12.06.2010 




sp. 


ECCUS 265 


JF901863 


Sivas 


30.05.2010 




tabidus (Fabricius, 1775) 


ECCUS 266 


JF901869 


Icel 


12.04.2009 




tabidus 


ECCUS 267 


JF901870 


Canakkale 


06.05.2010 




troglodyta (Fabricius, 1787) 


ECCUS 268 


JF901864 


Zonguldak 


03.05.2010 


Hartigia 




linearis (Schrank, 1781) 


ECCUS 270 


JF901896 


Ardahan 


07.06.2010 




linearis 


ECCUS 271 


JF901897 


Kirsehir 


03.06.2003 




linearis 


ECCUS 272 


JF901898 


Kirsehir 


03.06.2003 




nigra (M. Harris, 1779) 


ECCUS 273 


JF901899 


Konya 


17.05.2009 




sp. 


ECCUS 274 


JF901900 


Sivas 


17.05.2010 




sp. 


ECCUS 275 


JF901901 


Sivas 


13.05.2010 




xanthostoma (Eversmann, 1847) 


ECCUS 276 


JF901902 


Zonguldak 


03.05.2010 




xanthostoma 


ECCUS 277 


JF901903 


Zonguldak 


03.05.2010 


Syrista 




parreyssii (Spinola, 1 843) 


ECUUS 280 


JF901906 


Sivas 


26.05.2007 




parreyssii 


ECUUS 281 


JF901907 


Adana 


05.06.2003 


Characopygus 




sp. 


ECCUS 290 


JF901895 


Icel 


13.04.2009 


Pachycephus 




cruentatus (Eversmann, 1847) 


ECCUS 300 


JF901904 


Sivas 


06.06.2009 




smyrnensis J .P.E.F. Stein, 1876 


ECCUS 301 


JF901908 


Edirne 


07.05.2010 




smyrnensis 


ECCUS 302 


JF901909 


Edirne 


07.05.2010 




smyrnensis 


ECCUS 303 


JF901910 


Sivas 


11.06.2010 




sp. 


ECCUS 304 


JF901905 


Sivas 


12.06.2010 



and a 5 min final extension at 72°C. The purification and sequencing of amplification 
products were performed using a commercial sequencing company (Macrogen Ltd., 
Seoul, Korea.). Sequencing reactions were carried out in both directions using the 
same primers as in PCR reactions. The forward and reverse nucleotide sequences were 
assembled and edited by eye using the CodonCode Aligner v 3.5.6 (CodonCode Cor- 
poration) and aligned by using CLUSTAL W version 1.83 (Thompson et al. 1994), 
using the default parameters of the program. Finally, all the sequences obtained are 
deposited in GenBank (Table 1). 
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Data analysis 

Estimates of evolutionary divergence analyses were conducted in MEGA5 (Tamura et 
al. 2007) using the Kimura 2-parameter model (Kimura 1980) over sequence pairs be- 
tween genera. The rate variation among sites was modeled with a gamma distribution 
(shape parameter = 0.87). The presence of substitution saturation was determined with 
DAMBE version 4.5.18 (Xia and Xie 2001). The genetic distance versus the number 
of transitions and transversions at first, second and third codon position in all taxa was 
plotted to examine the saturation at a partial COI gene sequences. 

In order to investigate the phylogenetic relationship of Cephinae, phylogenetic 
trees were constructed using maximum parsimony (MP), maximum likelihood (ML) 
and Bayesian inference (BI) methods. Nucleotides were used as discrete and unor- 
dered characters. The best-fit model of DNA substitution and the parameter estimates 
used for tree constructions were chosen according to the Akaike Information Crite- 
rion (AIC) as implemented in Modeltest version 3.7 (Posada and Crandall 1998). 
The phylogenetic signal in the data partitions was estimated by maximum likelihood 
mapping method (Strimmer and von Haeseler 1997) using TREE-PUZZLE version 
5.2 program (Schmidt et al. 2002). MP phylogenies were estimated, with characters 
unordered and equally weighted, under the heuristic search algorithms 'simple' and 
'TBR' using PAUP version 4.0bl0 (Swofford 2002). Bootstrap estimates were cal- 
culated from 100 replicates under the above search options. This whole procedure 
was also applied to the data after removal of the third codon position. ML analy- 
ses (Felsenstein 1981) were conducted using RAxML-VI-HPC v. 4.0.0 (Stamatakis 
2006)f. The AIC results from Modeltest provided the GTR + I + G model as the 
best-fit for substitution model. BI analysis was performed with the software BEAST 
v. 1.5.2 (Drummond and Rambaut 2007). The analysis was run with four chains for 
5 x 10 7 generations, sampling from the chain every 5.000 generations. This generated 
an output of 10 4 trees. All analyses were performed assuming a Yule process of diver- 
sification. In order to confirm that the chains had achieved stationary, we evaluated 
"burn-in" plots by plotting log-likelihood scores and tree lengths against generation 
number using the software Tracer v. 1.5 (Drummond and Rambaut 2007). After 
determining convergence, we discarded all samples obtained during the first five mil- 
lions generations as "burn-in". The percentage of samples recovering any particular 
clade in a BI analysis represents that posterior probability of a clade (Huelsenbeck and 
Ronquist 2001). A majority rule consensus tree (Bayesian tree) was then calculated 
from the posterior distribution of trees, and the posterior probabilities calculated 
as the percentage of samples recovering any particular clade. The BI tree built with 
TREEANNOTATOR, discarding the initial 10% of samples as burn-in [Fig-TREE 
v. 1.3.1] (Rambaut 2008) was used to visualize the results. For the sake of a better 
presentation, branches representing individuals belonging to same species were col- 
lapsed if the species is recovered as monophyletic (Fig. 4). 
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Results 

Evaluation of the material collected after publication of Korkmaz et al. (20 1 0) and 
Calmasur and Ozbek (2010) revealed that there are two additional species of Cephus 
occur in Turkey. The examined material is presented below. 

Cephus parvus (Dovnar-Zapolskij, 1931) 

Material examined. Turkey: Sivas [39°42.71'N, 37°01.30'E] 1300 m, 26.05.2010, 
1$, 17.05.2010, \S. 

Distribution. Palearctic region. 

Cephus runcator Konow, 1896 

Material examined. Turkey: Edirne [40°39.32'N, 26°17.82'E] 50 m, 07.05.2010, 

6$,lc?. 

Distribution. Turkey, S. E. Europe. 

The complete alignment of the partial mitochondrial COI gene sequences from 
68 cephid specimens, including representatives of these two new records, resulted in 
a fragment containing 658 base pairs, among which 287 nucleotide positions were 
variable and 223 sites of which were parsimony-informative. The analyzed sequences 
correspond to a functional mitochondrial gene region because of the presence of sin- 
gular peaks in each chromatograph and absence of in-del and premature stop codons, 
and presence of the highest nucleotide substitutions at the third codon position (Avise 
1994). The percentages of nucleotide composition at each codon position are vari- 
able (Fig. 1). The mean frequency of COI sequences used in the analyses showed a 
bias of A + T (T 37.0%, C 15.2%, A 33.9% and G 16.0%), which is similar to other 
reported members of Hymenoptera (Jermiin and Crozier 1994, Dowton and Austin 
1995,1997, Leys et al. 2000, Danforth et al. 2003). The A + T content at the third, 
second and first codon positions are 90.7%, 59.8%, and 61.9%, respectively. The nu- 
cleotide G has lowest (1.0%) and the A highest content (52.8%) at the third codon 
positions. The distribution of polymorphic sites for all cephid species shows that the 
majority of substitutions are at synonymous sites. The vast majority of synonymous 
substitutions are also found at third codon positions with a rate of 87.88% for the 
Cephinae. The first and the second positions are relatively more conserved in compari- 
son with the third position. 

The numbers of base substitutions per site from averaging over all sequence pairs 
between genera are shown in Table 2. The least diverged genera appears to be Cham- 
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Avg 1 Pos 2 Pos 3 Pos 



Figure I . Percentage of nucleotide composition at each codon position. 

copygus and Pachycephus (p= 0.062) and, the most are Hartigia and Syristra (p= 0.161) 
also with highest standard error value of 0.017. 

All three codon positions in the partial COI gene were analyzed for saturation, 
achieved by plotting the number of observed substitutions against the model TN93 
genetic distance estimates. The scattergrams (Fig. 2a-c) showed that transitions and 
transversions for the first, second and third codons of the partial COI gene increased 
with the genetic distance, but considerable scattering was also observed. In addition, 
a similar plot of the third codon transition of the COI gene (Fig. 2d) suggested that 
saturation of transition occurred between certain pairs of the taxa, which may lead to 
higher levels of homoplasy (Kumar et al. 2001; Zhang et al. 2008). 



Table 2. Estimates of evolutionary divergence over sequence pairs between genera. The number of base 
substitutions per site from averaging over all sequence pairs between groups are shown. Standard error 
estimates are shown above the diagonal. 



Genera 


1 


2 


3 


4 


5 


6 


7 


8 


1. Calameuta 




0.010 


0.011 


0.010 


0.013 


0.012 


0.015 


0.025 


2. Tracheitis 


0.110 




0.010 


0.009 


0.013 


0.012 


0.014 


0.025 


3. Cephus 


0.108 


0.119 




0.007 


0.013 


0.010 


0.016 


0.026 


4. Chamcopygus 


0.078 


0.094 


0.069 




0.012 


0.007 


0.012 


0.028 


5. Hartigia 


0.136 


0.146 


0.143 


0.114 




0.013 


0.017 


0.029 


6. Pachycephus 


0.116 


0.125 


0.113 


0.062 


0.137 




0.013 


0.029 


7. Syrista 


0.124 


0.145 


0.156 


0.102 


0.161 


0.125 




0.030 


8. Outgrup 


0.248 


0.263 


0.256 


0.249 


0.292 


0.278 


0.279 
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Figure 2. Saturation plots of transversion and transition rates against JC69 distance at a first codon posi- 
tion b second codon position c third codon position, and d sum of data. 





Figure 3. Likelihood mapping analysis of the sequence alignments of COI gene present in the Cephinae. 
The regions at the corners of the triangles correspond to the three possible tree topologies for a quartet; 
the lateral regions to partly resolved trees and the central region to unresolved trees. The numbers indicate 
the percentage of quartets falling in each region. 
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Result of likelihood mapping is presented in Fig. 3. High dichotomic phyloge- 
netic signal was detected in the dataset. The percentage of the quartets suggesting a 
star- or network- like phylogeny is 9.9%, indicating that data are reliable for a di- 
chotomic phylogenetic analysis (Schmidt 2009). For ML analysis GTR+I+G models 
showed a significantly better fit than the other less complicated models for the COI 
dataset. Maximum likelihood analyses under the same model of evolution resulted in 
topologies with InL = — 5570.6831 in RAxML, which were very close to the BI tree. 
Bayesian inference under the GTR+I+G model resulted in a topology with mean InL 
= - 5347.963. Posterior probability values from the BI were congruent with ML boot- 
strap support. ML and BI analyses generated a tree with almost the same overall topol- 
ogy (Fig. 4). Equally weighted parsimony analysis of the 287 parsimony-informative 
characters produced 12 most parsimonious trees with a length of 1065 steps (Homo- 
plasy Index = 0.608, Retention Index = 0.392 and Consistency Index = 0.392). These 
equally parsimonious solutions were due to differences in terminal branches. However, 
the branching pattern of bootstrap tree was comb-like and recovered almost no origi- 
nal branches. Considering that this may be due to many synonymous changes in the 
third codon position, we run an analysis excluding the third codon position from the 
data. The analysis produced 60 equally parsimonious trees with a length of 200 and 
the bootstrap application was also resulted with no support to branching pattern of 
original trees. This may be partly attributed to the nature of data and relatively a short 
sequence not sufficient to detect phylogenetic signal under parsimony interference. 
Therefore, we do not present any MP trees here. 



Discussion 

Currently, the Cephinae is divided into three tribes based on morphology and feeding 
habits of larvae. The recovered mitochondrial gene trees substantially conflict with the 
current taxonomic arrangement, particularly the tribe level. Trees constructed under 
ML and BI methods supported monophyly of each genus except Pachycephus but failed 
to recover monophyly of any tribes. However, it should be noted that monophyly of 
most genera were supported by low posterior values (Fig. 4). This is probably due to 
the strongly biased nucleotide composition and the saturation at the third codon posi- 
tion (Fig. 2). The BI tree suggests that the most basal clade of Cephinae is the genus 
Cephus making the Cephini paraphyletic with respect to rest of Cephini and other 
tribes. Occurrence of Syrista within Pachycephini rather than Hartigiini makes both 
tribes polyphyletic and paraphyletic respectively (Fig. 4). Otherwise, Pachycephini and 
Hartigiini appear as sister groups. However, we do not propose a new classification as 
the present phylogeny is generated from a single gene fragment. 

Evolution of phytophagy has occurred many times in insects, and is often accom- 
panied by a significant increase in rates of speciation (Mitter et al. 1988). Phytopha- 
gous insects are also notable for their high degree of host-plant specialization; probably 
over 75% of species feed only on members of one plant family (Bernays and Chap- 
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Figure 4. Bayesian interface tree based on the mitochondrial COI gene sequences of the Cephinae. Host 
plants are indicated in parentheses. Numbers at nodes indicate the posterior values. 
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man 1994), and many insect species feed only on a single plant species (Scheffer and 
Wiegmann 2000). Syrista which is considered in the tribe Hartigiini, occurred within 
Pachycephini clade (see Fig. 4) and this placement is questionable as larvae of Syristra 
feed on Rosa. However, if this placement is considered to be true than it suggests a host 
shift event from Papaveraceae to Rosaceae. Occurrence of Cephus at most basal clades 
also suggests two later shifts from Poaceae to Rosaceae and Papaveraceae (Fig. 4). Con- 
sidering relationships among genera and species of Cephini inferred from the present 
phylogenetic hypothesis, several host shift events are also evident. However, we are re- 
luctant for further discussion until a more robust phylogeny become available derived 
from analyses of several gene sequences of both nuclear and mitochondrial genomes. 
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